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We propose dynamical schemes to engineer coherent states of a mechanical resonator coupled to 
an ancillary, superconducting flux qubit. The flux qubit, when repeatedly projected on to its ground 
state drives the mechanical resonator in to a coherent state in probabilistic, albeit heralded fashion. 
Assuming no operations on the state of the mechanical resonator during the protocol, coherent 
states are successfully generated only up to a certain value of the displacement parameter. This 
restriction can be overcome at the cost of a one-time operation on the initial state of the mechanical 
resonator. We discuss the possibility of experimental realization of the presented schemes. 
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I. INTRODUCTION 



The last decade has been witness to incredible progress 
in the fabrication of high quality factor, mechanical res- 
onators (MR) in the nano-scale characterized by funda- 
mental frequencies in the 100 MHZ to 1 GHZ range. Con- 
trol of such systems in the quantum limit has emerged 
as a very active field of research with implications for 
ultra-sensitive sensing technologies, quantum informa- 
tion processing as well as fundamental tests of quantum 
mechanics [l|. A key required feature for applications in 
the quantum limit is the ability to engineer well-defined 
quantum states of the MR. The properties of a ther- 
mal MR can be manipulated through controlled quantum 
mechanical coupling with electromagnetic radiation. In- 
deed, this forms a powerful avenue to cool the mechanical 
systems, coherently control their evolution and measure 
them. 

Recently, various techniques have been proposed to 
cool the MRs to their motional ground states, such as 
bang- bang cooling Q, single-shot state-swapping cool- 
ing via a superconductor [s'], side-band laser cooling 
[3|- [Tl|. electromagnetically-induced-transparency cool- 
ing [l2|, dynamic dissipative cooling [l^, and through 
repeated projective measurements on an auxiliary flux 
qubit puj ]. In this paper, we consider protocols for 
preparation of pure coherent states, instead of ground 
states, of MRs from initial general states (e.g., thermal 
states). Such coherent states can be used as inputs e.g. 
for schemes to generate superpositions of spatially sepa- 
rated quantum states (so-called Schrodinger cat states) 
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that can be used to study quantum behaviour, as well 
as the transition from quantum to classical behaviour, in 
macroscopically populated (phonon) states of MRs. We 
note here that a scheme for generating Schroe ding er cats 
in a mechanical system was discussed in Ref.|15| in the 
case of a MR capacitively coupled to a Cooper-pair box, 
which assumes the availability of an initial coherent state 
without providing for an explicit method of engineering 
it. 

In contrast with the conventional resolved sideband 
cooling scheme where the MR reaches the steady state, 
our previous work shows that the MR cooling can be re- 
alized by repeated projective measurements on the aux- 
iliary flux qubit [ij . Here, instead of ground state cool- 
ing of MR, we propose a cooling scheme to put the MR 
directly into a coherent state by repeated random mea- 
surements. 



II. THE DRIVING PROTOCOL: RANDOM 
PROJECTIVE MEASUREMENT 

We consider here a doubly clamped MR coupled to 
a gradiometer-type superconducting flux qubit psl [Tij 
(the schemes presented here formally work for coupling 
to any two level system). The simplest such system can 
be described by the effective Hamiltonian j20| : 



H = Hr, 



huj, 



hg{a 



(1) 



Here Hg is the free Hamiltonian of the flux qubit, the 
second term represents the free Hamiltonian H^iji of 
the MR with fundamental frequency a;„i, while the third 
term describes the interaction between these two systems 
Hint- In the Pauli spin- 1/2 representation, the flux qubit 
Hamiltonian in the basis of equal superpositions of per- 
sistent currents |0) ± |0) is described by a bias energy 
term ~ a^, that is set to zero by trapping a half- flux 
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quantum in the superconducting loop [ISl . |23| , and a tun- 
nelling term between the persistent current states, i.e. 



Hq ~ —cTz. We shall consider the weak coupling limit, 
i.e. g << ijJm- In this regime, the rotating wave ap- 
proximation applies and the system is described by the 
Jaynes-Cummings model: 



Hjc = Hq + Hmr - hg{a'' a + a ) 



(2) 



{ax ± iay)/2. 



wher 

In a previous paper by some of us [ij], based on the 
system described by the above JC Hamiltonian ([2|), a 
feasible scheme for the ground-state cooling of a MR 
was achieved via projective measurements on an aux- 
iliary flux qubit which interacts with it. For the pro- 
tocols to generate coherent states, we shall require the 
MR to be displaced in space by a. To this end, let 
us assume a flux qubit-MR model whose orginal Hamil- 
tonian is given through a displacement transformation 
D{a) = exp(aat — h.c.) to the JC Hamiltonian ([2]) (up 
to a constant term) 



H 



effective 



Hjc ~a^/2muJ^{a + a^) - gaa^. (3) 



Initially, we assume that the flux qubit -l-MR system 
is prepared in the product state \g){g\ (8 pm, where 
\g) = (|0) - \0))/V2 is the ground state of the free 
flux-qubit Hamiltonian, while the MR is in some general 
state with non-zero overlap with the coherent state we 
want to prepare. The basic constituent of the protocols, 
we present, for driving the MR to a coherent state is a 
quantum-Zeno type effect: the flux-qubit -f MR is allowed 
to undergo dynamic evolution under the Hamiltonian Q 
which is interrupted by a series of randomly timed pro- 
jected measurements on to the ground state of the flux 
qubit O — \g) {g\^I. We show below that after N consec- 
utive measurements on the qubit yielding outcomes \g), 
the state of the MR will approach to the coherent state 
|q) with the fidelity increasing to 1 with increasing N. 

Notice that the effective evolution operator pertaining 
to the MR encompassing a single measurement at time r 
can be disentangled to the form 



(5|e-^^"°1.g) - D{a)v('\r)D{a)\ 



(4) 



where from Eq.© vj^^r) = {g\ e-'"'"^^/''' \g) and fur- 
thermore D{a) = exp(Q;(a^ — a)) is the displacement 
operator [2l|. After one measurement 

Ou(^\g) (.g| ®p„)c/tOt = \g) (gl^plJl) (5) 

= D{a)V('Hr)D{aypmD{a)Vq^('Hr)D{a)^ (6) 

Thus after N consecutive measurements of the 
flux qubit ground state performed at times {r} — 
{ti, T2, . . . ttv}, the density matrix of the MR is 

, _ D{a)Vq^''\{r})p'i^Vq^^''\{r})DHa) 



where Vj^'^r) = Vg^'\n)Vg^'\r2) . . . vI'\tm) and pf = 
D^a)pmD{a). The success probability of this process is 

P;(7V) =Tr (i?(a)^W({r})p:;f y;W({r})i?t(«)). 

(8) 

We first analyze the properties of the above suc- 
cess probability. The Jaynes-Cummings evolution 
exp{—iHjcT/h), due to the symmetry associated with 
conservation of total number of excitations a^a -I- |e)(e|, 
where |e) ^ |0) -I- |0) is the excited state of the flux 
qubit, is easily decomposed into 2 dimensional sectors, 
apart from the the space of excitations which is 1- 
dimensional [13, H^l • The conditional evolution operator 
Vg^\T) therefore is diagonal in the phonon basis {\n)} of 
the MR (corresponding to Hmr). 



g 



Xn\n){n\ 



(9) 



with complex eigenvalues Aq e ^"^/^ and for n > 1, 
A„ = e-^("-i/2)<^'"^(cosl7„r + isinf7„rcos26l„) where 
ri„ — y%A^^cj^^JJ^y4T|r^2^ jg the Rabi frequency. The 
parameter 0„ measures the ratio of the effective inter- 
action energy scale to the deviation from resonance, ac- 
cording to the equation: tan20„ — 2gy/n/{A—ujm). Note 
that |Ao| = 1, is conserved due to effective decoupling of 
the zero phonon state during the conditional evolution, 
while for random times r, |A„| < 1 for higher phonon 
states n > 1. 

Hence after N consecutive measurements of the flux 
qubit ground state, the MR is driven into the state (0 
given explicitly as 



E 

n,fc^O 



A„(r)A^(7 



{n\p-i^\k)D{a)\n) {k\ {a) 



(10) 

where A„(t) = 11^=1 ^n{Tj). The success probability of 
realizing the conditional evolution (|5]) is thus: 



P;(7V)= Yl UryXl{n\p';^m4D{a)\n){k\D\a)) 

= E i^«(t)iV™^"^ (11) 

where the second line comes about from the cyclic prop- 
erty of the trace and unitarity of the displacement oper- 



ator D{a) and p' 



!ff(n) 



n\p^\n) is the diagonal matrix 



element of the effective input density matrix. 

On the other hand, the fidelity of the state of the MR 
PU)) with the coherent state \a) — D{a)\0) is given by: 



FUN) 



^ K{T)Xlir) 



D{ay pl,{N)D{am 
n|pf|fc)(0|n)(A:|0) 



n,k>0 



(12) 
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where in the last hne we have used the fact that 
|Ao(t)P = 1. This result shows that in order to obtain 
a non-zero fidelity with a coherent state a, the effective 
input density matrix must have a non-zero overlap with 
the phonon vacuum. 

We now turn to the asymptotic characteristics of 
our method. Due to the properties of the eigenvalues 
Xn{T) mentioned earlier, the products over different times 
|A„(t)| — > as the number of measurements — > oo for 
all n > and only |Ao(t)| 1. Thus asymptotically, 



p:{n) = p: 



cff (0) 



{a\p„ 



a 



(13) 



Using this, we see from Ea. (|T2| that the fidelity asymp- 
totically tends to 1, which shows that one effectively ob- 
tains the coherent state of the MR. Importantly, the re- 
sults (jl2ll3l) also imply, that by measuring the ground 
state of the ancillary system, in this case the flux qubit, 
many times, one basically remotely implements the co- 
herent state projective measurement on the initial state 
of the MR which is hard to implement directly. Indeed, 
Eq. (fT3]) shows that the asymptotic success probability 
of the protocol is given by the probability of projecting 
the initial state on the target, coherent state \a). 

The characteristics of the described process depend on 
the initial state of the MR. In what follows, we con- 
sider two, rather natural, possible, initializations of the 
MR and compare the performance of the above described 
protocol. 



III. THERMAL STATE OF THE MR 



D(^'")(a) = {i\ D{a) \n) reads 

= ^(i|D(a)at"|0) 
{^\ia^-ar\a) 



exp(-|a|^)V , 



X (-1)'^ 



min(n,?') 

E 

(i — m)! 



(14) 



In the last equality we had in mind that (a™ |i))^ = if 
m ^ i + 1. We also applied this property for coherent 
states 



{i - m\a) = exp(-|ap/2)- 



(15) 



After rearranging terms, this probability is written as 

N 

fM(^) = exp(-|an n \\l{r,)W) 

n i j — 1 



X n! i\ 



min(n.'i) 

E 



(16) 



m! {n — m) \ [i — m)! 



On the other hand, the fidelity ([T^ reads 



Consider the initial state of the MR to be the ther- 
mal state pm = exp(— /3fia;ma^a)/Tr[exp(— /3?iCLima^a)] = 
^j^p(n)|n)(n|, where the probability p(n) = r;,"/(l-f n)" 
with average phonon excitation at inverse temperature 
P given by the Planck formula n — l/[exp(/3?icj„i) -I- 1] 
and {\n)} is the phonon state basis of the undisplaced 
MR. This can be obained by allowing the MR to ther- 
malize, while the interaction with the flux qubit {Hint) 
as well as the "displacement" term for the MR, i.e., the 
second term in Eq. (jS]), is switched off. Then the inter- 
action as well as the displacement term are switched on 
at time t = Q followed subsequently by measurements of 
the qubit ground state. 

The success probability of N consecutive measurements 
([TT|) in this case becomes 



where p"(i) = '^^^ — "j^ is the phonon number prob- 
ability distribution in a coherent state a. 

We show below, that for this class of states, the MR can 
be driven into a coherent state only for coherent param- 
eters below some critical value Oc which depends on the 
temperature or equivalently the average phonon number 
of the thermal MR. Indeed, this stems from the decreas- 
ing overlap at a given temperature of the phonon number 
probability distribution p{i) of the thermal state, which is 
geometric with maximum value for the vacuum, and the 
phonon number distribution of the coherent state p"(i) 
which is Poissonian and centred at the average phonon 
excitation value . Thus to drive the MR into coherent 
states with higher average phonon numbers will require 
higher temperatures. 



, ^ ^ ^ IV. DISPLACED THERMAL STATES OF THE 

E ni^"(^^)l^'WI(*l^(")l")l MR 

For a given temperature, the restriction on the co- 
where the matrix elements of the displacement operator herent state parameter can be avoided, using a scheme 
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based on a different input state. Indeed, consider in- 
stead the thermal state corresponding to the displaced 
MR. Unlike in the previous scenario, we therefore con- 
sider a setup where the MR is first displaced, followed 
by a period of thermalization, after which the interac- 
tion is abruptly switched on followed by implementa- 
tion of the driving protocol. In this scenario, the input 
state is the thermal state of the displaced MR, p^'^ = 
exp [— /^ftwrn ((a^a — Q;(a -f a'f ))] /Z ~ D{a) pm D'' (a), 
where Z = Tr {exp [—/3?ia;,„ ((a^a — a(a -I- a^))] } and 
Pm is the thermal state of the undisplaced MR, as given 
in the previous section. 

Under this assumption, the new density operator ma- 
trix after N consecutive succesful measurements, is 



pM(iV) = |.g) {g\® 



j^(a)y/(r)p,„y;^(r)jt(a) 



(18) 



which is simply the displacement of the thermal state of 
the (undisplaced) MR. The success probability is simply 



PM(7V) 



(19) 



n>0 



which is independent of the coherent state parameter a 
and asymptotically (|13|) is (0|/5m|0), i.e. the occupation 
probability of the vacuum state of the bare MR at the 
given temperature. 



V. NUMERICAL RESULTS 

We now study quantitative characteristics of the two 
presented driving protocols numerically. We consider a 
realistic 27r x 100 MHz nano-mechanical resonator with 
quality factor Qm = 10^, coupled near- resonantly to a 
flux qubit non-resonantly with a tunnelling splitting A = 
l.lt^m- 

We first consider the scheme described in Sec. III. The 
MR is taken to be in equilibrium at the ambient temper- 
ature T = 20 mK, which corresponds to a mean phonon 
number of 3.69. The coupling constant is assumed 
g = 0.04a;m and the measurement times Tj are randomly 
selected. Figure [1] shows the evolution of the fidelity 
(|12p with increasing number of projective measurements 
for different values of a up to a < Ofmax — 5.1 above 
which the asymptotic success probability < 8 x 10^^. 
The fidelity approaches the asymptotic value of 1 for 
all a's within 30 measurement steps. In these cases, 
the MR is driven close to the coherent state \a) as can 
be seen also by considering the average phonon number 
(n) = 5n -|- |ap which is the sum of thermal-like and 
coherent-like contributions. The thermal-like contribu- 
tion 6n can be calculated by subtracting the coherent 
like contribution |ap from the average phonon number 
(n) in the calculated final state (fTO|). The thermal-like 
contribution after 30 measurements is n = 2 x 10~^. 




Coherent parameter 







N (measurements) 



FIG. 1: (Color online) Fidelity F^''>{N) for different values of 
a after A'' projective measurements using g — 0.04 LUm and the 
time interval within which a single projective measurement is 
performed is Tj = lO/ujm. 



In Fig[2l we consider the variation of the success prob- 
ability of obtaining N consecutive ground state measure- 
ments of the flux qubit. For a = 2 (red curve), this prob- 
ability changes very slowly after 20 measurements show- 
ing that the coherent state has practically been achieved. 
On the other hand, as we approach the value a = 5.1 the 
survival probability becomes very small, i.e. we practi- 
cally no longer can obtain the flux qubit in its ground 
state. Indeed for a = 5 (green curve) the survival proba- 
bility saturates to a very small value within ten measure- 
ments. 




10 15 20 

N (measurements) 



FIG. 2: (Color online) The success probability Pi^'(iV) for 
different values of a. Red triangles denote the case q = 2. 
Green triangles denote the case a = 5. The controllable pa- 
rameters are the same as in Fig. [T] 



From the results above, the asymptotic probability of 



5 



success (with number of measurements A'' — > oo) is seen 
to be a measure of obtaining a coherent state with fideUty 
1 . It decreases to zero with increase of the coherent dis- 
placement parameter a as shown in Fig. [31 



0.25r 




2 4 6 

Coherent parameter 



FIG. 3: (Color online) Asymptotic value of the success prob- 
ability depending on the displacement a that is sought to be 
obtained. 

We now consider the scheme described in Sec. IV. Fig- 
ured] presents the same physical quantities as Fig. [T] and 
Fig. [2] as well as the mean phonon number for a ther- 
mal state of the physically displaced MR. We consider a 
higher bath temperature T = 40 mK with mean phonon 
number 7.84. After 30 measurements, the thermal-like 
contribution is n ~ 10^"^, thus, the average phonon num- 



ber is (n) 



\a\ . It is noticeable that in this framework 



the coherent state is always reached after 30 measure- 
ments. As mentioned in the previous section, the char- 
acteristics of the protocol including the fidelity with the 
coherent state (see Eq. (fT8| ) is independent of the dis- 
placement parameter. 



when the flux qubit and MR are coupled both resonantly 
or off-rcsonantly. We showed that the driving process is 
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(c) 
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FIG. 4: (Color online) (a) Survival probability after A'' mea- 
surements, (b) Fidelity, (c) Themal-like mean phonon num- 
ber Sn{N) for the initial phonon number n = 7.84. The pa- 
rameters are T = 40 mK and r = 8/lj. 



VI. CONCLUSION 



akin to performing a projective measurement on the MR 
without directly acting on it. 



We would like to remark that though the results for the 
protocols of Sec. Ill and Sec. IV are obtained for the case 
of randomly-timed measurements, they are equally valid 
if one performs equal or fixed time measurements. Exact 
timing of measurements is however difhcult to implement 
experimentally and indeed unfavourable. As shown in 
[14| . randomly timed measurements result in faster dy- 
namics towards asymptotic states of the MR. 

In conclusion, we have demonstrated that quantum 
Zeno effect-type schemes consisting of repeated measure- 
ments on an auxiliary two level system can drive a nano 
mechanical resonator in to a coherent state. We propose 
two different ways to realize coherent state: the first one 
is efhcient for small values of a, while the second one, 
with a displaced initial state of the MR, is robust for 
any value of the coherent parameter. The schemes work 
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